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Abstract. Topological insulators have gap less edge/surface states with novel transport 
properties. Among these, there are two classes of perfectly conducting channels which are 
free from backscattering: the edge states of two-dimensional topological insulators and the one- 
dimensional states localized on dislocations of certain three-dimensional topological insulators. 
We show how these novel states affect thermoelectric properties of the systems and discuss 
possibilities to improve the thermoelectric figure of merit using these materials with perfectly 
conducting channels. 



1. Introduction 

Topological insulators (TI) [TJ [2J [31 01 [5] have attracted a lot of attention in the recent years. In 
these materials the excitations in the bulk are gapped just like in the normal insulators, but due 
to specific topological properties of the band structure these insulators possess gapless modes 
at the edges (on the surface). These modes are topologically protected and are stable against 
any perturbations that do not break the time-reversal symmetry. Various novel properties of 
TIs have been theoretically proposed, and some of them have been experimentally observed. 
The advent of this class of materials has brought us renewed interest in nonmagnetic insulators 
which have rarely been considered as an interesting subject for studying electronic properties. 
The transport properties of the edge/surface excitations in topological insulators are very distinct 
from those of the low-energy excitations of conventional metals. In two-dimensional (2D) 
topological insulators, the edge states are free from elastic scattering by nonmagnetic disorder. 
In three-dimensional (3D) topological insulators [6l[7], elastic scattering for the surface states is 
at least partly suppressed. It is also proposed that in certain classes of 3D topological insulators, 
the lattice dislocations have gapless states, which are robust against elastic backscattering by 
nonmagnetic disorder [8]. 

Thus, we have two cases of topologically protected one-dimensional (ID) gapless states in 
topological insulators; (i) the edge states in 2D topological insulators, and (ii) the states bound 
to dislocations in certain 3D topological insulators. These states have similar physical properties, 



that they are perfectly conducting, and they carry spin current. Hence these ID gap less states are 
expected to have large contributions to transport. Having in mind that some of the topological 
insulators are good thermoelectric materials, we studied the thermoelectric properties of these 
ID states, i.e. for edge states [9] and for states bound to dislocations [10]. In the present paper, 
we review the results of these articles and discuss them in detail. The paper is organized as 
follows. In Section 2 we study thermoelectric transport by the edge states in 2D topological 
insulators, following Section 3 is devoted for thermoelectric transport by the states bound 
to dislocations in 3D topological insulators [TO]- The results are summarized in Section 4. 

Thermoelectric properties of the system are characterized by the thermoelectric figure of 
merit ZT defined by 

ZT = I™, (1) 

K 

where a is the charge conductivity, S is the Seebeck coefficient, k is the thermal conductivity, 
and T is the temperature. The larger ZT means better efficiency for thermoelectric conversion. 
There are various approaches toward improving the ZT. In many materials, the main obstacle 
for better ZT is a large phonon heat transport. One approach to impove ZT is called phonon- 
glass electron-crystal (PGEC) [11], i.e. a good conductor for electrons but a poor conductor 
for phonons. The phonon conduction is suppressed by varying the lattice structure, and thus 
making the system glassy for phonons. Nevertheless the electron conduction should remain good 
for a good thermoelectric transport. Another approach for improving ZT is low-dimensionality 
|12] , This idea is to control the electronic band structure, so that the density of states has a 
peaked behaviour [13]. Despite these attempts the maximum ZT remains of order of unity, 
awaiting for a breakthrough. 



2. Thermoelectric transport by edge states of 2D topological insulators 

Because the 2D topological insulators have gapless edge states which are free from elastic 
scattering, they might have a possibility for novel and enhanced transport properties. For 
charge transport it has been measured experimentally in HgTe quantum well. In multiterminal 
measurement, the edge channels are shown to be perfectly conducting [161 Ej- While the edge 
states form perfectly conducting channels, the number of edge states is typically much smaller 
than that of the bulk, and the overall transport is mainly contributed by the bulk. In charge 
transport, we can neglect the bulk transport at lower temperatures, if the chemical potential is 
within the bulk gap. Nevertheless, because we are interested in the heat transport, we should 
consider finite temperatures, and bulk carriers cannot be neglected. Hence, in order to make the 
edge states more prominent compared with the bulk, we consider the system in narrow ribbon 
geometry. 

The charge and heat current in a solid is described by the following linear response, 



w J \ Li L 2 J\ -T~ l VT 



(2) 



where j is the current density, q is the electron charge, w is the heat current density, and fi is the 
chemical potential. Li are the transport coefficients, using which the various physical quantities 
are expressed as 

2 T a Li L\ \LqL 2 -L\ 

CT = qL , S = ——, 7T = — , K = K ph + K e , K e = — (3) 
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where K p h and K e are the thermal conductivity by the phonons and that by the electrons, 
respectively. 



2.1. Thermoelectric transport by edge states only 

There are three carriers which contribute to ZT: the edge electrons, bulk electrons, and phonons. 
We first consider an ideal situation of only edge states, while transport by the bulk electrons and 
phonons is neglected. The density of states is schematically described as in Fig. [D^a). We set 
the origin of the energy as the conduction band edge. We apply the Landauer-Biittiker formula 

L » = {fz2 UmiE - /*)" (- H) dE, (4) 

i 

where i is the channel index, 71(E) is the transmission probability for the i-th channel at energy 
E, and f(E) = \/(e( E ~* 1 '' knT + 1) is the Fermi distribution function [II]. In the present case, 
the transmission probability is equal to unity for the states with energy E at —A < E < 0, and 
zero otherwise, with A being the size of the bulk gap. We then obtain 

LI = ^ f° (E->,y(-%)dE = e e F e (v,»)(k B Ty, (5) 
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c e = -r, F e (v,fj,)= x v dx, (6) 
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where s is the cross section of the ribbon, £ is an effective length of the ribbon, A = A/k B T, 
and p, = fx/k B T. The superscript "e" denotes the edge transport. The effective length £ is taken 
as the inelastic scattering length of the edge states Anei- If we neglect the phonon transport, the 
figure of merit ZT is 

^ = T Ll ~ 2 (7) 

The resulting ZT is shown in Fig. Qlh). We can see that as the chemical potential becomes 
larger than zero and goes into the bulk band, the ZT increases rapidly. Compared with the ZT 
ever achieved in experiments which is around ZT ~ 2, the resulting ZT is extremely large. It is 
because ZT reflects the average energy per carrier in the unit of k^T; if the chemical potential 
is in the bulk band, the edge carriers have larger energy in average, giving large ZT. 

In this regime of fj,, the edge carriers are thermally excited holes. Therefore the carrier 
numbers are exponentially small, so are the transport coefficients L|. Hence, if we include bulk 
carriers or phonons, they will easily become dominant over the edge contributions, suppressing 
the large ZT. 

2.2. Competition between edge and bulk transport 

In addition to edge transport, we consider the bulk transport. The schematic figure is shown 
in Fig. [2j For simplicity, we set the bulk gap to be sufficiently larger than the temperature, 
and the chemical potential /i to be near the conduction band edge. Under these assumptions 
the bulk valence band is neglected. We also note that we are considering very narrow ribbons, 
for which the transverse motion is quantized due to the confinement. For ribbons as narrow as 
W = lOnm wide or so, we can assume that only the first subband is involved. 

Bulk transport is diffusive and calculated by the Boltzmann transport theory. For simplicity 
we put the lifetime r to be constant. 
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Figure 1. (a) Schematic density of states (DOS), (b) Figure of merit ZT as a function of 
chemical potential //, calculated from edge transport only. Bulk carriers and phonons are 
neglected. A is taken to be A = 0.15eV. 
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Figure 2. (a) Schematic picture for the two-dimensional ribbon, (b) Density of states including 
edge states and bulk states at the first subband. 

where e(k) = k 2 /2m, c is the number of carrier pockets, /i* is the mobility, m is the effective 
mass, and the superscript "b" denotes the bulk transport. The prefactors c e and Cb characterize 
the size of the transport for the edge and the bulk, respectively, and the total transport coefficient 
is given by their sum: 

L v = LI + L b u . (10) 

The results for the transport coefficients at T = 1.8K are shown in Fig. (a)-(d). We took 
Bi2Te3 as an example [28], and adopt the parameters as c = 6, £[ ne \ = 1/um, m* = 0.02m e . The 
cross section of the ribbon is set to be s = lOnmx 0.5 nm. The mobility is jx* = 2000cm 2 V- 1 s~ 1 , 
taken from the bulk mobility at 80K [18j . For the phonon thermal conductivity we assumed 
kl = 0.1Wm _1 K _1 estimated from the bulk values [191 [20] . 

In Fig. [3l we observe interesting interplay between the edge and bulk transport. It is 
reasonable that in the charge conductivity (Fig. (3]^a)) and the thermal conductivity (Fig. [3]^c)), 
when the chemical potential is in the bulk gap (fl <C 0), the transport is dominated by edge 
carriers. On the other hand, when the chemical potential is in the bulk band (jl S> 0), the 
transport is dominated by bulk carriers. The behavior of the Seebeck coefficient S (Fig. E]^b)) 
is interesting. First, the edge and bulk contributions have opposite signs. It is because the sign 
of the Seebeck coefficient is the same as the sign of the carriers involved. When the chemical 
potential is near the bulk band edge, the carriers in the bulk conduction band are electrons, 
whereas the carriers in the edge states are holes. Thus they have opposite contributions for 
the Seebeck coefficients. It is also noted that the Seebeck coefficient from the edge states is an 
increasing function of /i, while that from the bulk states is a decreasing function of [i. They are 
both asymptotically linear functions. 



As a whole, the figure of merit ZT (Fig. E^d)) has a maximum around /i = 0, the bulk 
band edge. We note that although the ZT from the edge is very large when /i is large, it is 
reduced drastically by the bulk carriers and phonons. It is because in this regime of /t, the 
carrier concentration of edge states is exponentially small, and it can be easily overcome by the 
bulk carriers and the phonons. 
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Figure 3. (a) Charge conductivity a, (b) Seebeck coefficient S, (c) thermal conductivity k, 
(d) thermoelectric figure of merit ZT. The horizontal axis represents p = fi/ksT, and p = 
corresponds to the band edge. 

We can evaluate the transport coefficients when p is far from the band edge, i.e. \p\ ^> 1. 
Here we also neglect the valence band and focus only on the conduction band and edge states, 
by setting A = A/(/cb^) — > oo. For the edge states, when /] > 1 (chemical potential is deep 
inside the bulk conduction band), the asymptotic form of F e (u,fj,) is 

F e (0,n) -e-P, F e (l,/t) ~-(i + £) e -P F c (2, M ) ~ (2 + 2/2 + p 2 )^ , (11) 

and when p <C — 1 (chemical potential is deep inside the bulk gap), 

F e (0,/i)~l-e* F e (l,M)~(£-l)e", F e (2, p) ~ y + (-/2 2 + 2/2 - 2)e^. (12) 

On the other hand, for the bulk states, when /] > 1 (chemical potential is deep inside the bulk 
conduction band), the asymptotic form of F^(v,fj,) is 

2 2 

F b (0,M)~VA> F b(l,M)~g^, F b (2,/x)~y^ (13) 
and when p <C — 1 (chemical potential is deep inside the bulk gap), 

F b (0, M )~^, F b (l,/i)~^(|-A), F b (2,/x)~^(^-3/2 + /2 2 ). (14) 



From these asymptotic forms, we evaluate the asymptotics for the transport properties, shown 
in Fig. [3] as broken lines. In particular, the ZT from the edge state only has the asymptotic 
form 



ZT e ~ (1 + A) 2 :/!» 1, 
Zr e ~A(i_£)2eP:/i«-l, 

and the ZT from the bulk state only has the asymptotic form 
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2.3. Universal behavior described by dimensionless parameters 
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Figure 4. ZTyq&x ci function of r and g. fj, is set to maximize ZT. 

We have calculated various thermoelectric properties using specific values of parameters. 
There is a wide range of choices for various parameters, such as sample geometry, material choice, 
temperature, and so forth. It is therefore important to develop a theory which is applicable for 
various choices of parameters. This can be achieved by introducing dimensionless parameters r 
and g defined by 
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The parameter r characterizes the ratio of edge transport to bulk transport, and g characterizes 
the phonon transport to bulk electronic transport. The figure of merit ZT is then expressed 
as a function of r, g and ft. We then maximize ZT as a function of p. In Fig. 0]we show the 
maximum value of ZT for each value of r and g. maximum. One can see from Fig. 2] that ZT 
has a minimum at r ~ 2.3. As we increase r, the edge states become more and more dominant, 
and bulk-to-edge crossover occurs. Because the edge and bulk tend to cancel each other, this 
crossover is accompanied by a minimum of ZT. In general the behavior of r and g as functions 
of the temperature T is complex. Nevertheless, the inelastic scattering length I is a rapidly 
decreasing function of T, as is seen from the case of quantum Hall effect with £- ine \ ~ y— 1-5 |21| . 



and we can estimate r ~ T~ 2 . Thus we conclude that by lowering temperature, edge states 
becomes dominant, and there occurs a bulk-to-edge crossover. The ZT thereby has a minimum 
and then increases again as a function of temperature. For lOnm-wide ribbon, the crossover is 
expected to be around 5K-10K. 

2-4- Discussion 

In the above estimate, the value of ZT is sensitive to the edge-state inelastic scattering length I. 
If I becomes longer, it is more favorable for edge-state transport since it is free from scattering up 
to a longer distance. It is difficult to qualitatively estimate t in topological insulators because it 
is dominated by electron-electron scattering. We can draw analogy with quantum Hall systems, 
where I is estimated to be l^m at IK [21J. On the other hand, in HgTe quantum well, known as 
a 2D topological insulator, perfectly conducting channels for the system size l//m are observed 
at 1.8K, meaning that I > lfim [16j . From these data we used I = l^m in our calculation. 

Because r is inversely proportional to the ribbon width, narrower ribbons are more favorable 
for edge-dominated transport, in agreement with our expectations. On the other hand, there 
is a minimum width for which the ribbon supports the perfectly conducting edge channels. It 
is determined by a penetration depth A of edge states. If the ribbon width is narrower than 
the penetration depth of edge states, there occurs a backscattering between the edge states on 
the opposite sides of the ribbon, and the edge state transport is suppressed. In HgTe quantum 
well the penetration depth is 6nm-50nm [221 124] . On the other hand, in some systems such as 
bismuth(lll) 1-bilayer films, A is predicted to be as short as a lattice constant [23} 124] . and the 
ribbon can be made very narrow with keeping the perfectly conducting edge channels. 

3. Thermoelectric transport by dislocations in 3D topological insulators 

Compared with 2D topological insulators where the edge states are free from scattering, the 
surface states of 3D topological insulators undergo some scattering. In a case of Dirac cone 
as seen in Bi2Se3, the perfect backscattering by 180 degrees is prohibited, but other scattering 
processes can happen. Therefore, the surface states do not show perfect conduction. 



(a) (b) 




Figure 5. (a) Dislocations in 3D topological insulator. The inset shows one such dislocation of 
the screw type, (b) Schematic picture of the band structure for the calculation. 

Nevertheless, it is theoretically predicted that in some of the 3D topological insulators there 
occurs a Kramers pair of gapless states on a lattice dislocation [8]. This state carries a pure 
spin current and is topologically protected. Thus these dislocations form perfectly conducting 
channels, as is similar to the edge states in 2D topological insulators. It is then expected 



that if we increase the number of dislocations, as shown in Fig. [5](a), the transport by these 
ID channels becomes dominant over the bulk transport and phonon transport. Increasing the 
number of dislocations in the crystal is also helpful for other reasons: the mean free path of 
phonons is reduced by disorder, and the bulk carriers localize below the mobility edge formed at 
the bottom of the conduction band. Thus more dislocations enhance the ID channel transport 
while at the same time suppress the bulk electron transport and phonon transport. 

The framework of the ZT calculation is similar to the 2D ribbon, and we do not provide it 
here. The details can be found in [ID]. All the calculations are performed for room temperature, 
T =300 K. We assume that the phonon thermal conductivity is reduced due to disorder from 
the bulk value of K p h = lWm _1 K _1 to K p h = 0.01Wm _1 K _1 for the average distances between 
dislocations ~ 3 nm. The gap between the valence and conduction bands is taken to be A = 0.15 
eV. We also set the mobility edge to be E m = 0.05 eV, measured from the bottom of the 
conduction band, see Fig. [5jb). The resulting ZT for the ID channel and the 3D bulk carriers 
are given in Fig. Eta) and (b), respectively. As a function of the concentration of dislocations the 
figure of merit ZT is calculated in Fig. E^c). The highest value ZT reaches when the chemical 
potential [a, which can be changed by an external gate, is near the mobility edge. Remarkably 
the maximum value of ZT is over 6 at room temperature. 

It is important to mention here that, if the impurities are non-magnetic, the transition to 
a bulk Anderson insulator should not destroy the ID gapless states, because the time reversal 
invariance is not broken. Note also that by increasing the density of dislocations and/or non- 
magnetic impurities one can reach even greater values of ZT. Nevertheless, these densities are 
less likely to be achieved and can as well result in the tunneling of electrons between the channels 
of neighboring dislocations. This tunneling could lead to the gap opening and localization of 
the protected ID states. This is the reason for our choice of an upper limit of the dislocation 
density of n ~ 10 17 m~ 2 which corresponds to an average dislocation spacing of around 3 nm. 

We note that not every topological insulator has dislocations carrying ID gapless protected 
modes. The condition for the ID gapless protected modes is expressed as B • = it (mod2-7r), 
where the vector ~M. V (so called time-reversal-invariant momentum) is a part of the Z2 topological 
number for the topological insulator considered and B is the Burgers vector of the dislocation 
[8]. Therefore, only the materials which Z2 topological numbers satisfy the above criterion 
can carry ID gapless modes through the dislocations. In this sense, the dislocations in Bi2Te3 
[25"1 [26j [271 EE] and Bi2Se3 do not carry ID gapless modes, but the dislocations in Bii_ x Sba; 
[291 [30] (0.07 < x < 0.22) possess ID gapless modes. 

4. Conclusion 

In the present paper, we consider two types of topologically protected ID gapless states in 2D 
and 3D topological insulators, as current-carrying paths for thermoelectric transport. In both 
cases, the ID modes alone give large thermoelectric figure of merit ZT, whereas the bulk states 
and phonons suppress the ZT. Thus in order to have good ZT, it is necessary to have high 
density of ID modes. It is accomplished by making narrow ribbons for 2D topological insulators, 
and by making dislocations densely distributed in 3D topological insulators. 
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